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\j  Abstract 

Using  the  Heat  Equation  to  formulate  the  notion  of  scale-space  filtering,  we  show 
that  the  evolution  property  of  level-crossings  in  scale-space  is  equivalent  to  the  max¬ 
imum  principle.  We  briefly  discuss  Altering  over  bounded  domains.  We  then  con¬ 
sider  the  completeness  of  the  representation  of  data  by  zero-crossings,  find  observe 
that  for  polynomial  data,  the  issue  is  solved  by  standard  results  in  algebraic 
geometry.  For  more  general  data,  we  argue  that  gradient  information  along  the 
zero-crossings  is  needed,  and  that  although  such  information  more  than  suffices,  the 
representation  is  still  not  stable.  We  give  a  simple  linear  procedure  for  reconstruc¬ 
tion  of  data  from  zero-crossings  and  gradient  data  along  zero-crossings  in  both  con¬ 
tinuous  and  discrete  scale-space  domains.  1/  «  '  >  , 


r 


1.  Scale-space  and  zero-crossings 

The  use  of  multiresolution  representations  is  an  important  idea  for  the  analysis 
of  signal  and  image  data.  Many  data  structures  have  been  studied,  including  Gaus¬ 
sian  pyramids,  difference-of-Gaussian  pyramids,  Laplacian  pyramids,  and  “scale- 
space”  formulations  [1,2,3].  The  latter  formulation,  to  be  described  briefly  below, 
can  be  used  as  a  continuous  model  of  the  other  formulations.  We  will  discuss  the 
representation  of  data  by  zero-crossings  in  scale-space,  and  consider  the  stability  of 
reconstruction  methods. 

The  natural  framework  for  the  analysis  of  scale-space  formulations  of  mul¬ 
tiresolution  representation  is  in  terms  of  the  heat  equation  [4,5].  Specifically,  let 
/( jc)  be  a  bounded  function  defined  for  xSlR".  (Arbitrary  dimensions  can  be  han¬ 
dled  with  little  additional  fuss  over  the  case  n  =  1;  we  will  later  comment  on  the 
case  of  a  bounded  domain  D  Q  1R"  ).  We  define  u(x,t)  to  be  a  bounded  solution  to 
the  heat  equation: 

du  . 

~r~  —  Am,  (Heat  Equation) 

ot 

u(x,0)  *  f(x). 

The  solution  is  given  by  convolution  against  the  fundamental  solution  to  the  Heat 
Equation,  which  for  the  domain  1R"  is  given  by 

u(x,t)  =  /  K(x-y,t)f{y)dy, 

IR" 


where 
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We  see  that  u(x,t )  is  obtained  by  blurring  fix)  by  increasingly  diffuse  Gaussians, 
parameterized  by  t>0,  with  standard  deviations  a  satisfying  2 <r  »  4t.  In  computer 
vision,  scale-space  sometimes  refers  to  the  (x,o)  variables  that  can  be  used  to 
reparameterize  the  domain  of  u.  We  retain  the  (x,t)  parameterization  to  keep  the 
linear  Heat  Equation  relation  for  the  function  u. 

Convolution  by  Gaussians  is  considered  special  for  many  reasons  [6,4,7].  We 
see  from  the  above  analysis  a  relationship  between  Gaussian  convolution,  the  Heat 
Equation,  and  the  Laplacian  operator.  Of  course,  Gaussian  convolution  enjoys 
other  properties;  for  example,  die  central  limit  theorem  implies  that  Gaussian  con¬ 
volution  is  easy  to  implement  by  an  iterative  procedure.  However,  we  also  see  the 
extent  of  the  similarity  of  difference  of  Gaussians  and  the  Laplacian  of  the  Gaus¬ 
sian;  namely,  since  K(x,t)  is  itself  a  solution  to  the  Heat  Equation, 

(AlO  (x,f)  -  -^( x,t )  -  lim(AT(x,f+T)  -  K(x,t))/r. 
at  t-0 

That  is,  the  difference  of  Gaussians  is  a  good  approximation  to  A  AT  as  the  separation 
between  the  spread  of  the  two  Gaussians  approaches  zero  (and  the  difference  is 
scaled) . 

Filtering  by  the  Laplacian  of  a  Gaussian  can  be  written  in  three  ways: 

A K*f~  JT*A/  =  A(AT*/). 

If  we  denote  the  result  by  v(x,/),  we  see  that 

(1)  v(x,r)  is  the  fix)  data  filtered  by  the  Laplacian  of  a  Guassian. 

(2)  v(x,r)  is  the  solution  to  the  Heat  Equation  with  initial  data  A/. 

(3)  v(x,f)  is  Au(x,0,  where  u  is  the  solution  to  the  Heat  Equation  with  initial  data 
fix). 

The  zero  set  of  v(x,r)  is  the  point  set  in  (x,f)  where  v  =  0.  The  set  might  be 
empty  (for  instance,  if  /is  subharmonic  or  superharmonic;  [8])  everything  (if  /  is 
harmonic;  [9]  ),  or  a  proper  subset  of  (x,f)  space.  In  the  latter  case,  zeros  can  be 
isolated  points,  lines,  and  surfaces  (but  never  regions).  We  distinguish  components 
of  the  zero  set  which  form  manifolds  of  codimension  one: 

Definition:  The  zero-crossings  of  v(x,f)  refers  to  the  point  set 

a{(x,/)|v(x,0  <0}  n  a{(x,r)|v(x,r)>0}.  ■ 

Zero-crossings  have  been  suggested  for  segmentation  of  imagery  by  edge  detection 
[10],  and  for  stereo  matching  and  motion  correspondence  between  pairs  of  images 
(e.g.,  [11]).  It  has  also  been  suggested  [7]  that  the  zero-crossings  are  a  nearly  com¬ 
plete  representation  of  A/.  Finally,  Witkin  [3]  observes  that  zero-crossings  in 
scale-space  evolve  as  r  increases,  and  are  never  created  at  some  nonzero  t  .  This 
property,  discussed  in  [7]  and  in  [6],  ensures  that  zero-crossing  surfaces  are  nested, 
one  within  another,  enclosing  regions  containing  the  face  {f  =  0},  or  forming  a  sheet 
meeting  the  face  { t  =  0}  and  extending  to  t  =  «.  The  property  can  be  given  a  pre¬ 
cise  statement: 

Evolution  property  of  zero-crossings:  Let  C  be  a  connected  component  of  the  set  of 
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zero-crossings  in  the  domain  {(x.OlxCR",  T1&&T2},  where  O^Ti  <  Ti.  Then 

c  n{(x,ob-ri}#0.  ■ 

In  the  remainder  of  this  paper,  we  wish  to  make  two  main  points.  First,  we 
establish  the  equivalence  of  the  ‘‘Evolution  property  for  zero-crossings”  and  the 
classical  maximum  principle  for  parabolic  partial  differential  equations,  thereby 
allowing  us  to  consider  bounded  domains  and  nonstationary  convolutions.  Second, 
we  consider  recon structibility  from  zero-crossings.  The  principle  result  is  that  if 
knowledge  of  the  location  of  the  zero-crossings  is  supplemented  with  gradient  infor¬ 
mation  of  v  at  only  those  points  in  the  zero-crossings,  then  there  is  a  simple  scheme 
for  reconstructing  some  of  the  data,  but  that  even  then  numerical  accuracy  of  the 
reconstruction  is  unstable. 

2.  The  Maximum  Principle' 

The  classical  maximum  principle  for  the  solution  to  the  parabolic  equation 
du/dt  -  A u  states  (see,  e.g.,  [12, 13, 14]): 

Maximum  Principle:  Let  D  C  1R"  be  open  and  bounded.  Suppose  u  is  a  solution  in 
T  **  {(x,r)  |  x€P,  0<t<T}  of  class  C2  which  is  continuous  in  the  closure  T.  Then 
u  assumes  its  maximum  at  some  point  (x,r)  for  which  either  x€dP  or  t  —  0.  ■ 

Next,  suppose  that  scale-space  construction  is  denoted  by  the  operator  v  *  Sg, 
which  is  to  say  that  the  scaled  function  v(x,t)  is  obtained  from  the  initial  data  g(x). 
In  the  previous  section,  we  defined  S  to  be  v(x,r)  *  K(-,t)  *  g,  where 
;(x)  *  A/(x),  but  we  can  imagine  more  general  operators.  In  any  case,  it  is  logical 
to  make  certain  assumptions  about  S,  although  all  we  will  require  is  that 

( 1)  If  g(x)  is  continuous,  then  Sg  is  continuous. 

(2)  S(-g)  »  -Sg  for  all  g,  and  if  v  «  Sg  and  g(x)  *  v(x,t),  then  v  =  Sg,  where 
v(x,r)  -  v(x,/+t). 

(3)  If  g(x)-0  as  |x  |-oo  and  v  *  Sg,  then  for  each  t,  v(x,r)-0  as  |x  |-oo. 

We  note  that  if  S  is  defined  by  convolution  with  Gaussians  as  in  Section  1,  then 
the  maximum  principle  holds  for  v  —  Sg  as  long  as  g  is  continuous;  further,  condi¬ 
tions  (1)  -  (3)  hold. 

Our  first  result  is: 

Proposition:  The  following  are  equivalent: 

(i)  The  maximum  principle  holds  for  solutions  v  =  Sg  using  continuous  initial  data 
g(x)  satisfying  g(x)-0  as  |x  |-oo. 

(ii)  The  evolution  property  holds  for  level-crossings  of  solutions  v  —  Sg  using  a 
scale-space  operator  S  satisfying  (1)  -  (3)  above  and  continuous  g(x)  satisfying 
g(x)-0  as  |x  l-oo. 

Proof:  We  first  show  that  the  maximum  principle  implies  the  evolution  property. 
For  if  the  evolution  property  fails  for  a  level-crossing  /  for  some  v(x,t),  then  by 
suitably  transforming  v,  g,  and  /,  we  can  assume  that  /a: 0,  and  that  there  is  a  solu¬ 
tion  v  *  Sg,  (with  g  — 0  as  |x|-oo)  with  a  component  C  of 
{(x,t)  |  0^fS7\  v(x,f)>/}  disjoint  from  the  plane  t  =  0.  Let  (xo,t 0)  be  a  relative 
maximum  in  C  .  Then  since  C  is  open,  there  is  a  bounded  cylinder  in  C  with 
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(xo,*o)  in  the  interior,  which  is  in  violation  of  the  maximum  principle. 

Conversely,  if  the  maximum  principle  fails,  then  for  some  cylinder  D  x[0,T] 
and  some  v(x,/)  given  by  v  -  Sg  with  g(x)-0  as  |x  |-«,  the  maximum  of  v  occurs 
either  in  the  interior  or  the  top  of  the  cylinder.  Either  way,  there  is  a  value  /  less 
than  the  maximum  but  greater  than  the  values  on  the  bottom  and  sides  of  the 
cylinder.  Thus  there  is  a  component  of  the  level-/  crossing  within  the  wedge  O&t&T 
which  lies  entirely  within  the  interior  of  the  cylinder,  and  thus  does  not  meet 
{i r  -  0}.  So  the  evolution  property  for  level-crossings  is  violated.  ■ 

We  illustrate  the  utility  of  the  proposition  with  three  observations.  First,  given 
the  equivalence  with  the  maximum  principle,  any  proof  of  the  evolution  property 
for  zero-crossings  that  does  not  either  use  the  maximum  principle  or  essentially 
redo  the  proof  of  the  maximum  principle  is  highly  suspect.  Since  the  maximum 
principle  is  slightly  delicate,  especially  in  the  absence  of  strong  regularity  assump¬ 
tions,  the  former  course  seems  more  appropriate. 

Second,  using  the  version  of  the  gradient  Hopf  maximum  principle  for  the  Heat 
Equation  [13],  it  is  not  hard  to  show  that  knowledge  of  the  zero-crossings  together 
with  gradient  information  along  with  the  zero-crossings  (or  even  just  one  zero¬ 
crossing  contour)  is  sufficient  to  determine  v(x,f)  uniquely  [IS].  In  section  3,  we 
give  a  more  constructive  discussion  of  this  point;  but  it  is  interesting  that  the  max¬ 
imum  principle  establishes  this  uniqueness. 

Finally,  we  see  that  any  scaling  method  obeying  the  maximum  principle  will 
yield  the  evolution  property  for  zero-crossings.  Under  fairly  severe  restrictions, 
this  leads  one  to  Gaussian  convolution  [6],  but  more  general  scaling  methods  are 
possible.  For  example,  blurring  by  a  parabolic  operator  of  the  form  du/dt  *  Lu, 
where  L  is  a  uniformly  elliptic  linear  second  order  differential  operator  with  non¬ 
constant  coefficients  will  certainly  still  give  a  maximum  principle.  In  fact,  L  can  be 
nonlinear  [16].  Moreover,  suppose  we  replace  1R"  with  a  bounded  domain 
P  C  1R",  and  insist  on  data  /(x)  with  compact  support  in  the  interor  of  D  .  We  may 
then  define  v  =*  Sg,  where  g  «  A/,  by  solving 

dvldt  a  Av  in  D  x(0,«), 
v(x,0)  =  g(x)  for  xiD 
v(x,t)  =  0  for  xidD  . 

This  scaling  is  not  given  by  convolution  against  a  Gaussian,  since  the  domain  is 
bounded,  but  nonetheless  obeys  a  maximum  principle,  and  gives  the  same  evolution 
property. 

3.  Completeness 

We  return  to  a  consideration  of  zero-crossings  of  data  v(x,r)  obtained  by  filter¬ 
ing  initial  data  / (x)  defined  for  x  €1R"  by  the  Laplacian  of  Gaussians,  AAT ( • , r ) .  The 
question  we  wish  to  address  is:  to  what  extent  do  the  zero-crossings  represent  /  (  r)0 
Clearly,  /(x)  can  at  best  be  reconstructed  to  within  an  arbitrary  additive  harmonic 
function  and  a  scalar  multiple.  However,  if  we  assume  that  /( r)-0  as  |x  |-oc.  then 
only  the  multiplicative  constant  is  of  concern. 
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Yuille  and  Poggio  [17]  make  the  observation  that  if  g(x)  and  hence  /(x)  is 
polynomial  in  x,  and  if  n— 1,  then  reconstruction  from  zero-crossings  is  theoretically 
possible.  They  also  refer  to  the  validity  of  the  observation  for  larger  n. 

We  note,  however,  that  when  g(x)  is  a  polynomial  in  x€R"  for  any  n,  then 
v(x,f)  is  a  polynomial  in  (x,/)€R"+1.  Accordingly,  the  zero-crossings  are  part  of 
the  analytic  varieties  of  the  polynomial  v  as  studied  in  algebraic  geometry.  It  is  well 
known  that  the  varieties  in  Cn  determine  the  complex  polynomial  defined  on  n  com¬ 
plex  variables.  It  is  not  as  commonly  used,  but  nonetheless  true,  that  an  n- 
dimensional  subportion  of  the  intersection  of  the  analytic  variety  with  R"*1  also 
determines  the  polynomial  [18].  Thus  the  case  of  polynomial  data  can  be  settled 
Ufith  algebraic  geometry. 

However,  since  the  determination  of  a  polynomial  by  its  varieties  is  essentially 
an  analytic  continuation  result,  stability  of  the  reconstruction  is  unlikely.  That  is, 
small  errors  in  measurement  of  the  zero-crossings  could  lead  to  arbitrarily  large 
errors  in  the  determination  of  g(x).  Put  differently,  there  can  be  widely  different 
initial  data  leading  to  nearly  identical  zero-crossing  data. 

Worse,  settling  the  case  for  polynomial  data  says  little  about  the  case  of  con¬ 
tinuous  initial  data.  Although  the  Stone-Weierstrass  theorem  says  that  a  continuous 
function  can  be  uniformly  approximated  by  a  polynomial  on  a  compact  set,  the 
zero-crossings  depend  on  the  initial  data  globally,  and  the  dependence  can’t  be  local¬ 
ized.  Further,  the  lack  of  stability  means  that  the  approximation  is  irrelevant.  The 
situation  is  similar  to  the  fact  that  a  polynomial  of  a  single  variable  with  all  real 
roots  is  determined  by  its  zeros,  but  that  given  all  the  zeros  of  a  continuous  func¬ 
tion,  one  knows  nothing  more  than  the  zeros. 

In  fact,  there  are  known  examples  of  pairs  of  functions  /i(x)  and  /2(x)  such 
that  the  corresponding  vi(x,r)  and  V2 (x,t)  have  identical  zero-crossings  at  all  levels 
of  resolution.  John  Daugman  supplies  the  example  (for  two  space  dimensions)  of 
/i(*i.*2)  =  sinxi,  and/2(x1,x2)  =  (sinxi)(2  +  cosx2). 

However,  if  the  zero-crossing  data  is  supplemented  with  knowledge  of  the  gra¬ 
dient  data  at  the  zero-crossings,  then  reconstruction  of  at  least  some  of  the  data 
g(x)  is  theoretically  possible  by  a  quite  easy  procedure,  given  below.  Details  of 
these  ideas  were  reported  earlier  in  an  unpublished  work  [15].  The  use  of  gradient 
data  for  the  representation  also  appears  in  [4],  but  the  gradient  data  there  is  not  lim¬ 
ited  to  the  zero-crossings.  The  use  of  gradient  data  along  zero-crossings  is  dis¬ 
cussed  in  [19].  Many  researchers  have  noted  from  a  casual  observation  of  zero- 
crossings  of  image  data  that  zero-crossings  with  large  gradient  magnitudes  are  of 
greater  significance  than  those  with  low  gradient  magnitudes. 

3.1.  Continuous  Case 

Specifically,  let  fl  be  a  bounded  connected  component  of 
{(x,r)|  *2:0,  v(x,f)#0},  and  denote  by  D  the  set  {x€R"|(x,  0)€n},  and  by  T  the 
zero-crossing  3nn{r>0}.  Let  t  be  a  value  such  that  T>sup{r |(x,r)  €fl}.  Next,  we 
set  g(x)  -  g(x)  for  xiD  ,  and  g(x)  =  0  elsewhere.  Finally,  let  b(x)  be  the  g<x) 
data  blurred  to  the  level  t: 
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My)  -  / K(y-x,r)g(x)dx. 

»• 

Using  Green’s  theorem,  it  is  easy  to  show: 

Proposition: 

My)  *  J*r(y-x,T-f)Vv(x,f)-iyi<T, 
r 

where  n  is  a  surface  normal  to  T  at  (x.r),  and  da  is  surface  area  measure.  ■ 

Thus  given  the  zero-crossing  T  and  Vv(x.r)  for  (x,f)€l\  then  the  blurred  data 
Mx)  can  be  constructed  by  a  simple  linear  process j  The  original  data  g(x)  can  be 
reconstructed  by  deblurring  the  b(x)  data  [20].  Deblurring  is,  of  course,  a  classic 
unstable  process.  The  situation  is  not  hopeless,  however,  since  g(x)  has  known 
compact  support,  which  might  be  used  to  advantage,  and  also  since  errors  that  occur 
are  predominantly  in  high  frequency  components,  which  might  not  be  as  essential  to 
visual  interpretability. 

The  lesson  of  this  section,  ultimately,  is  that  even  for  bounded  zero-crossings 
supplemented  with  gradient  data  along  the  zero-crossing,  reconstruction  is  still 
unstable.  We  defer  a  remark  on  relaxing  the  constraint  that  the  zero-crossing  be 
bounded  until  the  next  subsection,  where  we  consider  a  discrete  version  of  the  result 
of  this  section. 

3.2.  Discrete  Data 

For  simplicity,  we  treat  the  case  of  one  unbounded  space  dimension,  although 
the  results  extend  easily.  We  are  given  data  fif  i-  ■  •  •  ,-1,0,1,  •  •  •  ,  and  define 

1  ,  1  ,  ,  1  , 

gi  *  ■j/i  - 1  ~  y/i +  ~£Ji  * 1 ' 

We  define  the  filtered  data  vj,*  recursively: 

Vi.o  =  gi. 

1  ^  1  .  1 

yvit*+— V(  +  i,*. 

We  also  define  the  blurring  kernel 


Both  v  and  K  satisfy  a  discrete  version  of  the  Heat  Equation,  namely 

1  1  ,  1  , 

It  is  not  hard  to  prove  a  discrete  analogue  of  the  evolution  property  for  zero- 
crossings.  The  key,  as  one  might  suspect  from  Section  1,  is  a  discrete  version  of  the 
maximum  principle,  which  is  easy  to  establish. 

Let  (1  be  a  bounded  4-connected  collection  of  pixels  (t'./t)  with  a  nonempty  set 
D  =  {i|(i,  0)€ft}.  Let  T  be  an  upper  bound  r>max{<: \(i,k)  (fl),  and  define  b,  to 
be  the  data  g,,  HD  .  blurred  to  level  T: 
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Finally,  let 


bi  “  2  ki-J,T'8] • 

3(±l.0)fl  -  {</,Jfc)$n|(/±l,Jfc)«n}. 
3(0, 1)0-  {(i.Jfc) €n|(U  +  !)*«}, 


a(0,-i)n  -  {(i,it)€n|  *>o,  (i,jk-i)£n>. 

Then  simple  but  messy  algebra  allows  us  to  show 
Proposition: 

Vt  t  +  Vj+4  k 

€—1,1  (!,*)€  3(4, 0)n  z 


^i+€-./.r-k-i+*i-j,r-t-i 


[v/+t,k~vfi*] 


+  2  4vlik+i*A:<_j>T-k-i 

(1,*)  €3(0,i)fl 


2  4  VixKi-jj-k. 

(<.*)€  a(o,-i)fl 

To  reconstruct  data  by  the  above  equation,  choose  a  connected  component  of 
{(i.Jfc) | v (f.Jk) >0},  (or  respectively  <0  ).  If  the  component  extends  to  infinity  in 
either  coordinate,  truncate  the  domain  to  become  a  convenient  bounded  collection  of 
pixels,  and  denote  the  result  by  Cl.  We  store  the  sets  d(±i,o)ft,  d(o,±i)^,  and  D  as 
defined  earlier.  For  pixels  (i.Jfc)  in  3(1,0)^  (respectively  d(_i,o )Cl  ),  we  store  the 
information  and  vj+i,*,  (respectively  and  v,_it*).  For  pixels  (i.Jfc)  in  3(o,i)^. 
we  store  the  data  vif*+i  and  for  d(o,-i)n  pixels  we  store  vlt*.  Using  the  above 
equation,  we  choose  a  T  and  reconstruct  the  blurred  data  bj.  To  reconstruct  the 
data  gi  for  i€P ,  it  suffices  to  deblur  the  ft;  data  by  solving  for  g,  in  the  linear  equa¬ 
tions  defining  b(.  In  fact,  the  system  is  overdetermined,  although  still  poorly  condi¬ 
tioned,  especially  if  \D  |  or  T  is  large. 

In  order  to  make  the  computations  feasible,  it  is  necessary  to  modify  the  for¬ 
mulas  for  a  bounded  spatial  domain.  We  in  fact  solved  a  bounded  domain  problem, 
with  — JVsJisjV,  setting  vjt*  -  0  for  i  =  ±N.  The  blurring  kernel  K  is  changed  by 
this  modification,  but  the  proposition  carries  over  with  little  change. 


In  Figure  1,  we  show  a  1-D  signal  git  and  the  zero  crossing  separating  positive 
and  negative  regions  of  the  associated  v,  *.  Applying  the  above  procedure  to  the 
central  positive  component  fl,  we  obtain  the  reconstructed  /»,  data  shown  in  Figure 
2.  The  true  Jb,-  data  is  identical  to  essentially  machine  precision.  Finally,  using  the 
method  of  pseudoinverses  to  deblur  the  data  shown  in  Figure  2,  we  obtain  g,  for  i  in 
the  middle  range,  as  shown  in  Figure  3.  This  is  to  be  compared  with  the  true  initial 
data  in  Figure  la.  The  poor  correspondence  is  due  to  the  fact  that  the  deblurring 
problem  is  poorly  conditioned;  better  deblurring  results  are  obtained  if  the  amount 
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of  deblurring  is  very  small.  However,  this  requires  that  the  top  of  zero-crossing 
contour  enclosing  the  data  occurs  after  not  too  many  blurring  steps.  The  lesson 
learned  here  is  that  although  reconstruction  is  in  theory  possible,  practical  recon¬ 
struction  may  be  impossible. 
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Figure  la.  Initial  g(i) 


Figure  lb.  Zero-crossings 

Figure  1.  An  initiml  function  g,  and  the  zero  crossing  pixels  in  v,*,  where  v,  *  blurs  the  ini¬ 
tial  data  by  scaling  in  k. 


Figure  2.  The  reconstructed  function  fr,  using  the  formula  from  Proposition  2.  The  data 
represents  the  result  of  blurring  the  g,  data  restricted  to  the  central  positive  interval,  zero  ex¬ 
tended  elsewhere,  to  a  level  k  above  the  top  of  the  zero  crossing  for  v(  t  shown  in  Figure  1 
above.  The  reconstruction  uses  only  information  about  v,  *  along  the  zero-crossing,  and  is 
nearly  exact  to  machine  precision. 
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Figure  3.  An  attempt  at  reconstructing  g,  using  the  data  b,  from  Figure  2.  The  attempt  fails 
because  the  matrix  equation  Kg  =  b  relating  the  g  data  to  the  b  data  is  poorly  conditioned, 
even  though  there  are  many  more  b  values  than  g  values.  Thus  the  small  errors  due  to 
round-off  in  representing  the  bt  data  are  magnified  when  reconstructing  g,.  The  singular 
value  decomposition  software  in  Creve  Moler’s  "matlab”  package  was  used. 


4.  Comments 

Zero-crossings  of  scale-space  filtered  data  seems  like  an  unlikely  form  of 
representation  of  data.  The  results  presented  here  suggest  that  even  when  supple¬ 
mented  with  gradient  data  along  the  zero-crossings,  the  representation  is  still 
unlikely.  However,  the  instability  of  the  representation  does  not  completely  deny 
its  utility,  since  it  might  happen  that  the  classes  of  functions  mapping  to  similar 
representations  share  properties  essential  for  interpretation.  A  required  step  in  the 
validation  of  the  utility  of  a  representation  is  an  analysis  of  the  invariant  properties 
of  signals  that  yield  similar  representations.  A  necessary  condition  is  that  attempted 
reconstructions  differ  from  originals  in  unessential  ways,  (from  the  standpoint  of 
interpretation).  The  methods  outlined  here  should  prove  useful  in  verifying  or 
disproving  this  necessary  aspect  of  establishing  a  viable  representation. 
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